PINNING OF INTERFACES IN RANDOM MEDIA 



N. DIRR, P. W. DONDL, AND M. SCHEUTZOW 

Abstract. For a model for the propagation of a curvature sen- 
sitive interface in a time independent random medium, as well as 
for a linearized version which is commonly referred to as Quenched 
Edwards- Wilkinson equation, we prove existence of a stationary 
positive supersolution at non-vanishing applied load. This leads to 
the emergence of a hysteresis that does not vanish for slow loading, 
even though the local evolution law is viscous (in particular, the 
velocity of the interface in the model is linear in the driving force) . 



1. Introduction 

Problems of interface evolution in heterogeneous media arise in a 
large number of physical models. Common to such models is a regular- 
izing operator, for example line tension, and the competition between 
an external applied driving force F and a force field f(x,y) describing 
the inhomogeneities. Assuming a viscous law for the relation between 
the driving force and the velocity of the interface, an important ques- 
tion is whether rate independent hysteresis can emerge in an average 
sense from the interaction between the heterogeneous force field and 
the regularizing operator. 

In this article, we consider a model for the evolution of an interface 
driven by its mean curvature through a random field of obstacles. Let 
neN,n>l. Let (fl, F, P) be a probability space, to G Vt. We model 
the interface as the graph (x, u(x, t, u)) of a function u : W 1 x K x Q — y 
K. moving through a field f(x,y,u) of (soft) random obstacles and a 
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constant driving force F. More precisely, we consider the PDE 
(1.1) 

_ , . r ._ , ,. / Vu(x,t,u) \ 

o t u{x, t,cu) — y 1 + \Vu{x,t, co)\ div — 

\n^l + \Vu(x,t,uj)\ 2 J 

+ ^l + \Vu{x,t,u)\ 2 (/(ar, u(x, t, u)) + F) 

=: J 1 + \Vu(x, t, <jj)\ 2 (k(u(x, t, u)) + f(x,u(x,t,co)) + F) , 

(1.2) 

u(x, 0, uj) = 0. 

The first term on the right hand side of equation ( 11. ip is the mean cur- 
vature operator for a surface that is given as the graph of the function u. 
The second term is the driving force, split up into the non-homogeneous 
random part / and the external constant loading F. The random field / 
will be specified in Section |2] in Condition 12. II and 12. 21 Basically we will 
assume that the nonhomogeneity consists of randomly distributed indi- 
vidual obstacles of a fixed smooth shape and possibly random strength. 
By n(u(x,t,u})) we denote the mean curvature operator for the graph 
of a function u(-,t,u) evaluated at x. 

Equation (11.11) is motivated in the following way (see also |CDL10j ): 
A very basic model for an interface (phase boundary, dislocation line 
in its slip plane etc) moving through an array of random obstacles (e.g. 
impurities, other dislocation lines) in an over-damped limit (inertial 
effects are neglected) is the gradient flow of the area functional plus 
a random bulk term. Consider thus a bounded set U C M n+1 and a 
smooth hypersurface V that is the boundary of the set A-p C U and 
define the energy 

(1.3) E(T):=n 1 (T)+ [ (f(X,u) + F)dX. 

J A T 

Here, H l denotes the 1-dimensional Hausdorff measure. The first inner 
variation (i.e., deforming the interface with the flow of a smooth vector 
field) yields the mean curvature k and a bulk term. The viscous gradi- 
ent flow with respect to the energy ( II .3p is thus given by the evolution 
equation 

(1.4) v n (x) = K(x) + f(x,w) + F, xeT, 

for the normal velocity v n of the interface. Now it is also possible to 
extend the notion of an evolving interface to U = M n+1 . The model is 
called quenched, because the random field does not explicitly depend on 
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time. For forced mean curvature flow and applications, in particular 
in the case of periodic forcing, we refer to |CLS09[ ICB04[ IDKY08] . 
Since we are only interested in constructing a supersolution for (11.41) 
when the initial surface is flat, it is sufficient to consider the mean 
curvature flow for an interface that is the graph of a function u, i.e, 
r = {(x,y) : y = u(x),x G R n }. 

If the gradient of u is sufficiently small, the evolution by forced mean 
curvature flow (MCF) for the graph can be approximated heuristically 
by a semi-linear parabolic PDE of the form 

(1.5) d t u(x, t, cu) = Au(x, t, uj) + f(x, u(x, t, cu), cu) + F rail™, 

(1.6) u{x,0,u) = 0. 

These kinds of problems have found considerable interest in the physics 
community, see e.g. |Kar97t IKle04t IBN04] . They are often referred to 
as the Quenched Edwards- Wilkinson (QEW) model. 

The goal of this article is to construct, for some F > 0, a stationary 
supersolution v to (jl.ip and to (II. 5p satisfying v > 0. In this article, 
we consider the natural case where, due to the randomness of the ob- 
stacle field, there exist arbitrarily large simply connected domains with 
positive combined driving force (i.e., where f + F > 0), in other words, 
large areas without obstacle. This makes a purely analytical approach, 
as employed in |D Y06t [DKY08] . insufficient. 

In order to illustrate the difficulty, consider in 1+1 dimensions a peri- 
odic array of obstacles with a forcing F > 0. Now remove each obstacle 
independently with rate p < 1 and increase the obstacle strength by 
S(p) in order to keep the expected value of the obstacle strength equal 
to the periodic one. The resulting random configuration (periodic with 
random "holes") may have a solution which is unbounded as t — > oo: 
For any h G Z and I G N, there exists almost surely az6Z such that 
the ball B((z,h),l) of radius I centered at (z, I) is free of obstacles. 
For sufficiently large I, the solution in this ball will grow to a height 
that is larger than h + 1. For an appropriate choice of obstacles and 
parameters, this perturbation can grow like a kink-antikink pair in a 
reaction diffusion equation (e.g. Allen-Calm) with forcing, until the 
entire curve has moved up at least one unit. Then the process repeats 
itself at a "hole" at height level h + 1. 

If such a supersolution exists, by the comparison principle for the 
mean curvature flow and for parabolic equations, an evolving solution 
u with any initial condition below v will always remain below v — the 
interface is pinned. Such pinning of an interface leads to a hysteresis 
that does not vanish for slow loading in the physical system. To see this, 
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consider a loading cycle starting with F = 0, increasing at first. The 
interface remains pinned until the driving force reaches a critical value 
(see Section [3] for a brief discussion of depinning). Above the critical 
force the material transforms (switching polarization, for example). 
Upon reversal of the driving force, the same phenomenon occurs^. One 
can see that the plotting the transformed region vs. the driving force 
will show a hysteresis loop that does not vanish even for slowly varying 
driving force F. 

In the following section, we construct such a supersolution under 
suitable assumptions on the random obstacle field. Section [3] provides 
some outlook. 



We first pose two conditions that will fix the structure of the non- 
homogeneous term / in equations (11 .11) and (11. 5p . This random nonlin- 
earity / is constructed in the following way: We consider an obstacle 
function <p £ C°°(IR ri x R) with the properties 

Condition 2.1 (Obstacle shape). There exist ri,r with r x > y/nr > 
0, so that 

i) (j) < 0, 4>(x,y) = for \ \{x,y)\\ > n, 

ii) 4>{x,y) < -1 for IKx,?/)!!^ < r . 

This fixes a 'shape' for the individual obstacles. Here, ||-|| denotes 
the Euclidean norm on M n+1 , II- II denotes the maximum- norm. The 

" II I I oo 

heterogeneity / is now given as the sum over individual obstacles with 
centers {(xi(co), y^cu))}^, and strength fiiu) > 0, i.e., 



We now pose a condition on the distribution of obstacles. The con- 
dition basically states that there is a uniform lower bound for finding 
an obstacle of some strength (also bounded from below) in a box of 
volume 1, independent of its shape or position, and independently for 
pairwise disjoint boxes. 



There is a difficulty in modeling this behavior. In a physical situation, the obsta- 
cles (non-transforming inclusions, for example) will always obstruct the evolution 
of an interface moving in any direction. A reasonable way to express this would be 
to consider the driving force f(x,u(x,t,uj),uj) ■ sign(<9 f w). Such an additional non- 
linearity in the equation would make the analysis unnecessarily complicated — we 
thus restrict ourselves to treating the transformation (dtu > 0) and the back- 
transformation (dtu < 0) separately. 



2. Construction of a supersolution 
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Condition 2.2 (Obstacle distribution). The random distribution of 
obstacle sites {(x k ,yk)}keN C M n x [r^oo) and strength {fk}keN C 
[0,oo) satisfy 

i) (xk,Vk) are distributed according to ann + 1-dimensional Pois- 

son process on R n x [7% 00) with intensity A > 0. 
H) fk ^re iid strictly positive random variables which are indepen- 
dent of {x k ,y k }. 

Note that there are no obstacles crossing the line {y = 0}, so at 
t = the interface with initial condition ( 11. 2ft only sees the external 
driving force. For a small time, the velocity of the interface is thus 
uniformly positive for F > 0. The comparison principle ensures that 
we thus have dtu > for all times. To see this, assume that for a 
smooth solution to (II. 5p there exists a first time to > when there 
exists xo G M. n : d t u(xo,to) = 0. Differentiating (11. 5p with respect to 
time yields 

d u u(x ,t ) = Au t (x ,t )+f u (xo,u(x ,t ))d t u(xo,t ) = Au t {x ,t ) > 0. 

Noting that, for a smooth solution of (jl.ip . the spatial gradient of d t u 
also vanishes at (xo,to), one can obtain the same non-negativity result 
for mean curvature flow. 

As mentioned in the introduction, the main difficulty in this work 
stems from the fact that there exist, with positive probability, arbitrar- 
ily large areas with arbitrarily small obstacles. We thus rescale space 
into boxes calling them "open sites" if they contain a suitably large 
obstacle. Theorems 12.41 and 12.151 below then depend crucially on the 
existence of an infinite cluster of open sites in Z™ +1 that is the graph 
of a Lipschitz function, for site percolation with P{Site open} > p c 
independently for all sites. It will be clear that the Lipschitz condi- 
tion is essential for the construction of a supersolution. This result is 
proved in [DDG + 10] , we repeat it here for the reader's convenience. 
See also [GH10J for an improved estimate on the critical percolation 
threshold. 

Theorem 2.3 (Dirr-Dondl-Grimmett-Holroyd-Scheutzow). Let n > 1 

and p G (0, 1). We designate z G Z n+1 open with probability p, and 

otherwise closed, with different sites receiving independent states. The 

corresponding probability measure on the sample space Q = {0, 

is denoted by P p . We write for the 1-norm on Z n+1 . The following 

holds: 

For any n > 1, if p > 1 — (2n + 2) -2 =: p c then there exists a.s. a 
(random) function L : Z n — » N with the following properties, 
i) For each x G Z n , the site (x,L(x)) G Z" +1 is open. 
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ii) For any x,y G Z n with \\x — y\\i = 1 we have \L(x) — L(y)\ < 1. 
Hi) For any isometry of Z" the functions L and L o have the 

same laws, and the random field (L(x) : x G II 1 ) is ergodic 

under each translation of HI 1 , 
iv) There exists A = A(p, d) < oo such that 



where v = (2n + 2) (1 — p) < I. 

We first show the result asserting the existence of a stationary pos- 
itive supersolution to the semilinear equation (jl.5p . since the calcula- 
tions are somewhat simpler. 

Theorem 2.4 (Existence of a pinned solution for QEW). If Condi- 
tions \2.1\ and \2. 6 A are satisfied, then there exists F* > and a non- 
negativev: l n xO — > [0, oo) so thatO > Av (x, u)+f(x, v (x, u), u)+F* 
a.s., i.e., any solution to (11. 5p with initial condition (11.61) and F < F* 
gets pinned. 

The proof consists of a piecewise construction of the supersolution, 
so it is first necessary to give some estimates on the components that 
will be used. We denote by B r the open unit ball of radius r around 0. 

Definition 2.5 (Local solution). Given r out > r in > ; F in > 0, and 

Font < 0, let v in be the unique solution of Av in = F in on B r . n C M. n , 
v in = on dB r . n . Let v out be the unique solution of Av out = F out on 
B rout \ B r . n C M n , with boundary conditions v out = on dB r . n and 
v ■ Vf out = on dB rout . 
We define 



Proposition 2.6. The function v\ oca i defined above satisfies 

i) ^locai is radially strictly increasing on B rout , 

ii) the graph ofv\ oc& \ restricted to B r . n is contained in the set B r . n x 



P p (L(0) > k)< Av k , k>0. 





otherwise. 



L 2n ' in' U J 



Hi) given f > 0, i 



(2.7) 
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and z/max{r in , f^-r^} < r , then v\ oca \ satisfies, in the sense of 
distributions (and in the sense of viscosity solutions), 

(2.8) > Aw local + f </>(., v local (-) +r ) + F on B rout 

for all F < min{-F out , / - F in }. 

Proof. The individual assertions are proved by a simple calculation. 

i) Follows immediately from the maximum principle. 

ii) The function v- m is nothing but a parabola, namely v- m (x) = 
4^ \x\ 2 — l^-r? . The assertion can be read off this form. 

Hi) From ii) and from the assumption on <fi in 12.11 one can see 
that for x G B r . m , we have <f)(x, V\ oca \(x) + ^o) < —1- The prop- 
erty (12. 8p for each individual piece of fi oca i can then be seen 
directly from the definition of t>i OC ai- The assertion follows by 
noting that equation (12.71) implies that the first derivative jumps 
down going radially outward across dB r; E This implies that 
"^locai is the pointwise minimum of two supersolution, thus a 
supersolution itself. The Laplacian is then a negative measure. 

□ 

Definition 2.7 (Rescaling). Given I > 2r\, d > 0, h > 0, and k = 

(k u k 2 ,...,k n ) e Z n , j G N, let 

V Qk ■= X n i=1 \ k S + d) + n, h(l + d) + l- n], 

H) Qk ■■= X i=l [h(l + d), h(l + d) + I], 
m ) Q ■= \J k Q k and D := R n \ Q, 

iv) Qk,j = Qk x [(j - l)h + n,jh + n]. 

Here /\ denotes the cartesian product of the n intervals following. 

Remark 2.8. The sets Q and D split M n into cubes, each separated 
by a distance d. The reduced cubes Qk are smaller by the length 2r% in 
every dimension, so that an obstacle with center in Qk fits completely 
inside Q k - The sets Qkj are extended in the n + 1-st direction by a 
height h. 

Proposition 2.9 (Percolating obstacles). Given h > 0, fix 1(h) > 
and f > s.t. 

l-exp{-A|A|-P{/ 1 >/}}>p c 

for\A\ = (l-2 ri ) n h, i.e., 1(h) = C (p c ,A,P{/i > f})h- 1 ^ n +2r 1 . Then 
there exists a random function L : Z n — > N with Lipschitz constant 1, 

2 The term — r- m is the radial derivative of v in and — F ont + r °^ 1 ) is the 

radial derivative of v ut at dB Tin 
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such that, a.s., for all k G Z n , there exists i G N such that (xi,yi) G 

Qk,L(k) and ft > /" 

For each k G Z n we select one obstacle index i e N with the above 
property and collect these obstacle indices in the set I. 

Proof. This is a direct consequence from Theorem 12.31 Indeed, consid- 
ering a cuboid Qkj open if it contains an obstacle of strength greater 
or equal /, we find that a cuboid is open with probability greater than 

Pc □ 

Definition 2.10 (Flat supersolution). We define the flat supersolution 
u flat : lR n ->■ E as Wfl a t(x) := min ie/ Vi ocal (x - X;). 

Proposition 2.11. Fix h > 0, f > 0, and 1(h) as in Proposition ^. 9[ 
Let r out = ^Jn{l{h) + | — r{) and assume that r out , r in , F m and F out 
satisfy the conditions in Proposition \2.b\ Then t>fl at satisfies, a.s., in 
the sense of distributions (and in the sense of viscosity solutions) 

> Au flat (x) + M x ~ x ^ ^flat(^) + r ) + F on R" 
for all F < min{-F out , / - F in }. 

Proof. Since ffl at is a minimum over shifted copies of the function twa, 
which is a supersolution where it is not equal +oo as proved in Propo- 
sition [521 it is enough to show that VQ at (x) < +oo for all x G M. n . This 
is, however, true by the choice of r ou ^ = y/n(l(h) + | - ri) with the 
property that the union over all k G Z n of closed balls of this radius 
with centers anywhere in cover all of M n . □ 

Remark 2.12. Since the function fi oca i is strictly increasing on B rout , 
the minimization process assigns each obstacle center Xi its Voronoi 
cell. On the Voronoi cell associated with Xi, the function fi OC ai centered 
at Xi attains the minimum. 

Proposition 2.13 (Gluing function). Fix h > 0, d > 0, I > 0. Let 

L : Z n — > R be a function with the property that if x,y G Z™ with 
\\ x ~ 2/11 1 = 1 we have \L(x) — L(y)\ < 2h. Then there exists C\ > 0, 
depending only on the dimension n, such that there exists a smooth 
function v g i ue : R n — > R such that 

i) for all k G Z n ; v g \ ue (x) = L{k) if x G Q k , 

H) SUpp Vfglue C D, 
III) IP^giuell^ < Ci£, 
iv) HVUgiuell^ < C X \. 
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Proof. It suffices to take a piecewise constant function that changes 
values on the center hyperplanes of the set D and apply a standard 
mollifier of size d/2. □ 

We now collected all the components to construct the supersolution. 



Proof of Theorem \2.4\ First, fix / and the function 1(h) as in Proposi- 
tion [279J Then fix < F in < | and r in such that max{r in , ^rf n } < r . 

We now need to find d > and h > 0, such that, with r out chosen as 
in Proposition 12.111 



(2.9) F in r in > |F out | ( -r in + ^ 

-nn + (v 7 ^ [1(h) + \ - rijj 

= l^outl (-r in + -1_ ^ (a,/^ 1/n +n + ^ 
and, at the same time, 

(2.10) |F out | >2C 1 A>2||A %ue || oo . 

Putting the two together, and using the fact that F in and r in are now 
fixed, it is sufficient to choose d and h so that 

(2.11) c ' >c Ti( h ~ 1/n + t + ri 

Since one has C± (h~ l l n + § + n) n < 2 n C + £ (f + n) n ) , one 
can see that there is such a choice. Now we fix F out = —2Ci-jk. 

Now we choose, according to Proposition 12.9} the index set / of 
relevant obstacles. From Proposition 12.131 and the Lipschitz condition 
on the percolating cluster of selected boxes from Proposition 12.91 we 
know there exists a function t> glue whose derivative is only supported 
on the set D and the property v g \ ue (xk) = yu + r for all k G I, and 

l|A%ue|L<^A. 

Choosing < F* < min{— |} One can now see that the function 

(2.12) V = U flat + Wglue 

satisfies 
(2.13) 

> Av{x) + ^2f(f){x-x i ,v{x)-y l ) + F* > Av(x) + f(x,v(x),u)+F*. 

□ 
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Remark 2.14. For a slightly different model, if the pinning sites are 
centered on a regular lattice, i.e., f(x, y, lo) = J2iez n jez+1/2 /ijC^O'K 3 ' - 
i,y—j), there is a lower bound for h — one can not make the boxes more 
shallow than the lattice spacing. This leads to the fact that there might 
not exist a d satisfying the estimate (12. lip . For n = 1, one can still 
find the supersolution in the described way, since the scaling of the 
gradient of v out with the distance d works favorably. In particular, the 
construction works for the model used in |CDL10j . 

For n > 2, the construction only works for / sufficiently large. De- 
pending on distribution of fi, such a choice for / might not be possible. 

We now turn towards the construction of a supersolution for the 
mean curvature flow. The theorem itself remains unchanged. 

Theorem 2.15 (Existence of a pinned solution for MCF). If Con- 
ditions \2.1\ and \2.2\ are satisfied, then there exists F* > and a non- 
negative w: R n xQ — > [0, 00) so that a.s., > k(w(x, u))+f(x, w(x, 
F* a.s., i.e., any solution to ( II. ip with initial condition ( II. 2p and 
F < F* gets pinned. 

The idea of the proof is to construct a local solution and a gluing 
function for the mean curvature operator and then provide estimates 
akin to Propositions 12.61 and 12.131 for these functions. The rest of the 
proof, modulo an estimate for the behavior of the nonlinear mean cur- 
vature operator when adding the local solution and the gluing function, 
can remain unchanged. 

Definition 2.16 (Local solution for MCF). As in Definition \2. 51 fix 
r in G (0,r ), r out > 0, F in > 0, and F out < 0, but now making sure that 
r- m < Fm and \F out \ is sufficiently small so that r i r °"r nn l". ^ n > |-F out |. 

'out I' out 'm; 

We construct the local solution from rotationally symmetric surfaces of 
constant mean curvature, so called Delauney- Surfaces |Del41] . 

Let w- m : B r . n -»■ R be given as w- m (x) = -\J F? n - \x\ 2 + yj Ff n - rf n . 
The radially symmetric function w out : B r . out \ B r . n — > R is defined by 
an elliptic integral as 




where C = L 




dp 
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We define 

w- m on B r . n , 

w out on B rout \ B r 

lim r ^ rout w out (r) on dB Tou 



W\oca\ 



out 5 



oo otherwise. 



Proposition 2.17. We have u>i OC ai is finite on B rout . Furthermore, it 
holds that 

(2.14) > re(w loca i) + /</>(-, ^iocai(-) + r ) + F on B rout 

in the sense of viscosity solutions, if 
i)0>F in -f + F, 
a) o > F ou t + F, 
m) F m - y/F? n - r? < r . 

Proof. The first statement is clear by inspection, since under the condi- 
tions on r in , F in , r out , and F out the functions w in and w OVLt remain finite. 
The second statement holds due to i) and ii) on the inside of the sphere 
B r . n (Condition Hi) confines the graph of Wi n to the set where < — 1) 
and on the inside of the annulus B r , \ B r . . Condition iv ) ensures that 

' out \ ' m / 

the derivative of w\ oca \ jumps only downwards going radially across the 
boundary from the sphere to the annulus, so that the mean curvature 
of wi oca i at the boundary is negative in the viscosity sense. □ 

Proof of Theorem \2.15l In order to employ the construction from the 
proof of Theorem 12.41 we first need to make sure that the scaling of 

\d r w ont (r)\ r=r . n \ = 



(rout-n n ) 2 " 2 
('•out-(rW-r in )") 2 F 2 u 



is suitable. Consider thus F out = — c^— and note that \d r w out (r) 
is then given by 

i 

g(r out ,c) : = 



(rout-r in ) 2 "- 2 r^ _ 
(r " ut -(r out -nn)") 2 C 2 r 2 r 2 



One can see by a simple calculation that for c small enough there 
exists C*2 > so that g(r out ,c) < C^c. This, however, implies that in 
the correct regime the Delaunay- Surface from the construction of w ovAi 
admits the same scaling properties with respect to r out as the function 

^out- 
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The proof of Theorem 12.151 can now be completed by the same con- 
struction as for the semilinear equation. First, fix the supersolution 
inside the obstacles. This determines the maximal outgoing radial 
derivative d r Wi n (r)\ r=r . n := G and the radius of the inner sphere r in . 
Consider then the function w flat constructed analogously to t>fl at above. 
It is necessary to satisfy (after setting r out = y/n(l(h) + | — ri)) 



g(r out , c) <G 

as well as 

(2.15) |F out | = c- s > 2Ci — . 



The scaling property discussed above ensures that this is possible. 

It remains to show that adding the function t> g i ue from Proposi- 
tion 12.131 does not destroy the property of negative mean curvature. 

Define v(v) := y 1 + |Vit| 2 . We have, after collecting terms from ex- 
panding the divergence in the mean curvature operator, 

A(u7 out + Uglue) 



^(Wout + Vglue) 



v{w mt + fgi uc ) 

^gluc ^gluc 11 



+ 



«(WoutJ 



3 

ue y 



K w out + fglue) 

1 1 



+ Am 



K^out + %ue) f(tWout) 
(D 2 W out ■ Vw out , VWout) 



+ 



K^out + %luc) 3 ^(^out) 3 
(D 2 

^glue ^gluG ), Vw out + V^gl 

uc 7 

K™out + % 
2(£> 2 W out ■ Vw out , V%ue) 

K W out + %ue) 3 
(£> 2 W out ■ V% ue , VVglue) 
v(w out + %ue) 3 
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K(w out ) + 0(||At; glue ||J + OiWAw^tW^ \\Vv g 1,2 



+ O(|| J D 2 «;out|| 0O i|V«; out ||^||V %ue || 2 J 

+ ^(|| D \ luc 1 1 ^ ( 1 1 Vw out + V %ue \\ 2 J 

+ 0(||DVu t || 00 (||V %U e|| 00 ||V Wout || 00 +||V^ lue ||L)- 

Here, (•, •) denotes the scalar product in MP and D 2 v ■ Vw denotes the 
matrix- vector product of the matrix of second derivatives of v : M n — > W 
applied to the gradient vector of w : M. n — > M. Note that HVugiueH^ = 
0(5), H^giucIL = 0(±) (Proposition M), and HVu^l^ = O(c). 
One can see that also | \D 2 u OVLt \ 1^ = 0(c), since for small curvature 
the gradient term dominates when calculating the second derivatives 
of the function t> out . All error terms can be made small with respect 

r n-l 

to \K(w ont )\ = \F ont \ = c (^ (Co/t -i% + |_ ri) y'M by noting that one can, 

instead of (I2.15p . for a given C > 0, fix h and d so that — F out > Ck. 
The rest of the proof then follows that of Theorem 12.41 □ 



3. Conclusions 

We have shown that, for our models of interface evolution in ran- 
dom media, a finite critical force is required to propagate the interface 
through the body. Many questions in this area, however, remain open. 
It was shown in |CDL10j . that for a model with obstacles on lattice 
site^l with random exponentially distributed strength for n = 1, no 
more stationary solution can exist if the forcing exceeds a critical value. 
The question whether interfaces in this case move with a finite speed 
of propagation is still open and currently under investigation (this is 
of course trivial for uniformly bounded obstacle strength when also 
avoiding overlap of obstacles). These two results together would show 
that there is a transition from a viscous kinetic relation (i.e., veloc- 
ity = force) in the microscopic model turns (after a time-rescaling) into 
a rate independent model for the macroscopic behavior of the system. 
Such rate independent kinetics are commonly assumed in macroscopic 
models of phase transformations or plasticity. This article provides a 
step into deriving this assumption from microscopic viscous kinetics. 



3 As pointed out in Remark 12.141 our construction of a supersolution for suffi- 
ciently small external driving force also works in this 'lattice case' for n = 1. 
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